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Abstract: We study /3-deformed matrix models with Penner type potentials, which cor- 
respond to M = 2 SU{2) supersymmetric gauge theories with Np = 2, 3, and 4 flavors. 
We compute explicitly the genus one corrections to the free energy of the matrix model 
and show that they match the corresponding results obtained from the Nekrasov partition 
function. 
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1 Introduction 

The past several years have seen much progress in the study of = 2 super confer mal field 
theories (SCFTs). It was shown that the compactification of M5-branes on a Riemann 
surface with punctures gives rise to a class of = 2 SCFTs [1]. Matter hypermultiplets 
of these SCFTs are related to the puctures on the Riemann surface. Elementary S-duality 
transformations are related to different sewings of the same Riemann surface. 

Moreover, it was proposed by Alday, Gaiotto and Tachikawa (ACT) [2] that the 
Nekrasov partition function of a class of N=2 SU{2) quiver SCFTs is identified with a 
chiral half of the correlation function of Liouville field theory on the corresponding Rie- 
mann surfaces. Specifically, the perturbative part of the partition function corresponds 
to the Liouville three-point function of primary operators and the instanton part of the 
partition function coincides with the conformal block generated by the Virasoro algebra. 

In [3], the ACT relation was generalized to A/" = 2 SCFTs with gauge group SU{Nc), 
relating them to the ^at^-i conformal Toda theory. [4] Also, the instanton partition func- 
tion of pure super Yang-Mills, which is the non-conformal limit of SCFTs obtained by 
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sending the mass of the matter hypermultiplets to infinity, was identified with the norm of 
a coherent state, called Whittaker state, in Liouville field theory. [5], [6] Numerous works 
have been done investigating various aspects of this surprising relation. 

In particular, Dijkgraaf and Vafa (DV) [7] observed that the classical spectral curve of 
a Penner type matrix model can be interpreted as the Seiberg-Witten curve [8], [9] of A/" = 2 
SOFT. They geometrically engineered the gauge theory with a suitable local Calabi-Yau 
geometry and employed the large A'" duality of the B-model topological string to describe 
the gauge theory with the matrix model. Given the CFT description of the matrix model, 
where the number of insertion of screening operators in Toda theory is the rank N of the 
matrices, [10], [11], it was suggested that the matrix model bridges between the gauge 
theory and the Liouville/Toda CFT. Matrix models corresponding to gauge theories with 
less number of flavors are also proposed in [13] . These matrix models in the context of the 
AGT relation have attracted much interest. [16]-[31] 

According to the DV proposal, the coupling gg and the Q background parameter e = 
ei + £2 = gs{-\fP — are of the same order in Qs- Thus, we double expand the free 

energy in gg and e, fixing the ratio Q = e/ gg or equivalently the central charge of Liouville 
theory. [12] In other words, we consider the free energy as 



In this paper, we will be computing the genus one part Fq^2 of the free energy in the /3- 
deformed matrix model. It has been shown in [14] that the genus zero free energy of the 
matrix model agrees with that of the SU{2) gauge theory. This evidence was extended to 
the half genus case of the /3-deformed matrix model in [15]. The genus one correction in 
the ordinary matrix model, i.e. f3 = 1, was studied in [19]. 

In section 2, we give a quick review of the /3-deformed matrix model with logarithmic 
potentials and derive the spectral curves from the loop equation. In the following three 
sections, we explicitly calculate the genus one part of the free energy of the matrix models 
for Np = 2, 3, and 4. In section 6, we summarize and discuss the results. In Appendix 
A, we present the detailed computation of integrals used in obtaining the free energy for 
Np = 3 and 4 cases. In appendix B, the Nekrasov instanton partition function is reviewed 
and the free energy is calculated from them. 

2 /3-deformed Matrix Model and Spectral Curve 
2.1 Penner type Matrix Model 

The partition function of /3-deformed matrix model is given by 



k,l>0 



— Fo,o + e-fo,i + e^-fo,2 + giFifi + • • • • 



(1.1) 




(2.1) 
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where A = Y\j^j{^i — ^j) is the Vandermonde determinant. The ^1 background parameters 
are related to the matrix model parameters through 

ei=5sV^, ^2 = -^. (2.2) 
The potentials for Np = 2, 3,4 cases are given by 

Viz)N^=2 = (2//3 + e) log ^ + As (^z + , (2.3) 

V{z)np=3 = {2i23 + e) log z + 2mi log(z - 1) - ^ , (2.4) 

V{z) Np=4. = (2mo + e) log z + 2mi log(2: - 1) + 2m2 log(z - q) , (2.5) 

where the mass parameters uiq, mi, 1712 with additional related to the four anti- 

fundamental hypermultiplet masses of the gauge theory by 

m = mi + ruoo, Ai2="2i-moo, fi3 = m2+mQ, ^^ = m2-mQ (2.6) 

The A2 and A3 are dimensionful parameters and correspond to the dynamical scales of the 
gauge theory, whereas q is dimensionless and is identified with the exponential of the UV 
coupling of the SCFT. The potential for Np = 3 case is derived from Np = 4 potential 
by taking //4 — )• 00 with fixed A3 = g/U4. The potential for Np = 2 is obtained from the 
Np = 3 potential by sending ^2 ^ 00 while = A3^2 is fixed. The neutrality condition 
in Liouville theory states that 

fii+fi3 + 9sVPN = 0. (2.7) 
The free energy of the matrix model is defined to be 

F = g^,logZ=-eie2logZ, (2.8) 
and it can be expanded in gg and e = ei + £2, 

F= Y^gfe'F.^i. (2.9) 

k,l>Q 

In ordinary matrix model, where /? = 1, the parameter e vanishes so only those terms with 
/ = survive. In this paper, we compute the leading e^-correction Fo,2 to the free energy. 

2.2 Spectral Curves from Loop Equation 

We define n-point connected resolvent by 

The loop equation, which is obtained from the variation of the partition function (2.1) by 
5\i = -j^—^ for small a, is given by [15], [32], 

glW[z, z) + W{zf + e W\z) + V\z)W{z) - f{z) = , (2.11) 
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where the function f{z) is 

/M-..v^E( ''''1:r'' )- <^-'^' 

We expand the resolvents and the potential as follows. 

W{zi, ... ,Zn)= 9se^Wk,lizi, ■■■ ,Zn) 

k,l>0 

Viz) = Vo{z) + eViiz) 
Then, the loop equation (2.11) is expanded to 

gy : Wloiz) + V^iz)Wo,oiz) - f{z) = 

gy : 2yo,o(^)^^o,i(^) + <oW + ^i'W^o,o(^) = (2.13) 
gy : 2yo,»{z)Wo,2{z) + Wl^{z) + Wl,^^{z) + Vi{z)Wo,i{z) = 

up to e^, where we have defined 

yo,o(^) = Wo,,{z) + \v:,{z) = \^V^{zY+^f{z), (2.14) 

which is the leading order term in the double expansion of the spectral curve, 

y{z) = W{z) + \v'{z) = ^ g^VkA^) . (2.15) 

k,l>0 

It is easy to see from the expansion of the loop equation (2.11) that Wifi{z) and 1^1^1(2;) 
vanish. 

The filling fraction in the matrix model is identified with the Coulomb branch param- 
eter a of the gauge theory and is given by the integral of the one form y{z)dz along the 
A-cycle of the spectral curve, 



fc,/>0 



where we have double expanded the vev. To compute the period integral, we first express 
yo,iiz) and 2/0,2(2) in terms of yofi{z) = yo{z) and the potential. From (2.13) and (2.14), 
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we obtain 



1 

Wo 



W^,,{z) + W^4z) + V({z)Wo,iiz) 

= -i(^o. + .;a-^<-^W'). (2.18) 
Using the expression (2.17) of 2/0,1(2:), we write yo,2(-z) as 



Thus, the subleading terms of the vector multiplet vev read 

ao. = ^<fyl^^dz, (2.20) 
27ri/ 2yo ' ^ ^ 

1 //y? (^o' + W)' , V{' + 2V{' 



= — N^ - ^ " ' i^^^ + ' ' ' \dz, (2.21) 
where we have dropped from ao,i the shift caused by the total derivative of logyo- [15] 
3 Np = 2 Model 

We first consider the matrix model which is dual to the SU (2) gauge theory with Np = 2. 
The action is given by 

V{z) = Vo{z)+eVi{z) = 2/i3logz + A + +elogz, (3.1) 

where we have omitted the subscript from A2 in (2.3). The function f{z) in (2.12) is 
evaluated to be 

m = 7 + 3 ^ (3.2) 

with 

ci = 5. (j2 (-^^77^ ^ ^) ) " 9sV^NA = -(/.I + ^3)A , (3.3) 
C2 = gsVP^(j2j^^ ' (3-4) 



- 5 - 



where the equation of motion V'{Xj)) = and (2.7) have been used in computing ci. 
The planar spectral curve (2.14) becomes 

yo(^)' = ^^, (3.5) 
where the quartic polynomial P^{z) is given by 

P^{z) = z +-^U3 + C2- — Iz - — z + (3.6) 

From now on, we will take fii = fi^ = m for simplicity. Then we have 

P,{z) =z^- —z^ + —z^-—z + l (3.7) 

with A = m^ + C2 — A^/2. The leading order vector multiplet vev ao,o = and the half 
genus contribution ao,i in (2.20) are written as 



1 / A WPJz) , , , 

1 / 1 / \\ dz l9ao 

We are going to express the genus one part 09,2 of the vev in terms of the derivatives 
of Co, which are given by 

dz 



dao 


1 


OA 


~ 27d 


d'^ao 


1 


dm? 


~ 2TTi 


d^ao 


1 


dmdA 


~ 2TTi 


d^ao 


1 


dA^ 


~ 2TTi 



(3.10) 



The ao,2 in (2-21) consists of three terms and we will consider the three terms in turn. The 
first term is 

2^/|3^^=2^/a(i6^-2^ + 7^J'^^- ^'-''^ 
The first integral on the right hand side can be written as 



2m J 16A p5/2 2m J 16A \ 3 dz"^ JPl 3 p3/2 



' ^ +^ dz, (3.15) 



27ri J 6A \^s/Pi 4p3/2 



-6- 



where we have integrated by parts in the last equahty. The exphcit form of P^iz) in (3.7) 
satisfies 

z^P^ = 3zPi + ^{z^ + z)- -^z^ . (3.16) 

Given the above relation and using the derivatives of oq in (3.10), (3.12), (3.13), we find 
(3.15) to be 

1 / ^2 Idao ^ I r z P' , m d^ap , Ad^ap 



2TTi J 16A p5/2 6 dA 2m J 8 A p3/2 4 dmdA ?, dA^ ' 
Thus, we obtain the first term (3.14) of 09,2 as 



= - A — -I - A * —dz 

27ri 7 8y3 6 OA A dmdA 3 dA^ 2Tri J 8A p3/2 



12 OA 4 am9^ 3 dA^ 
where we have used in the last step 



5 dap m d'^ap A d'^ap 



' ^ - %dz = -L / J^dz = 2A^ . (3.19) 



2m J p3/2 27ri / ^/Pl dA 
The second and third terms of ao,2 in (2.21) are given by 

1 / W + 1 / 1 , 2 , -,^2 



o dz = d) — -(z + 1)" — ^rn^dz 

2-Ki J 32y|] 2'KiJAK^ ' p3/2 

- 

8 (9m2 ' 

1 / W2 + 2Vl' , 1 / / 1 



27rz / 8yo 27rz / V 4Ay 

1 dap 

~ ~4 a4 ■ 

Plugging (3.18), (3.20), and (3.21) into (2.21), we find 



(3.20) 



(3.21) 



^Idap m d'^ap A d'^ap I d'^ap 

~ QdA^ A dmdA ^ 3 + 8 5m2 ' ^'^ > 

The planar vev (3.8) can be computed by doing elliptic integral and expanding hyper- 
geometric function, which was done in [14]. 

r-r f rn? ,0 {A^ - Qm^ A + Ibm'^) . fyiZm^ A^ - lAm'^ A + 2lm^) ^ 

ap = A 1 -tA^ - ^ ^A* -T, ^A'' 

" V 4^2 64yl4 256^6 



15(yl'^ - 2%m^A^ + 294^,4^2 - 924771^^ + 1001 



m" 



16384y48 



-A^ + 0(Ai°) . (3.23) 
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The half genus part (3.9) is given by [15] 

m 2 3m - 5m2) ^ 5 {sA^m - 28Am^ + 63m^) g 

105 (m (A^ - 21 A^m^ + 99Am^ - 143m6)) « 

^ — ^ TTF^ +■■■ ■ 3.24 

4096^15/2 ^ ' 

The genus one part (3.22) is 

A + m^ ^ A^ - 30^m^ - 35m^ 4 5 (3^^ - 49A^m^ + 147W + 231m6) 
«0'2 - - 16^5/2 A + ^^^^^§75 A i024yli3/2 

35 (^^ - 126^3^2 _^ 792^2^4 _ 358^^6 _ 2l45m8) „ 

+ — ^ TTTo ^A^ + • • • . 3.25) 

16384yli7/2 ^ ' 

We invert the equation a = + eao,i + £^00,2 to find A in terms of a: 

, n m? .9 — Qa?"rin? + 5m^ . ^ ba^m?' — lAa'^m'^ + 9m^ . = 

^ = a'' H ^A^ H A* H A'' 

' 2a2 32a6 64aio 

5a^ - 252a^m2 + 1638a^m'^ - 2860a'^m^ + U69m^ . . 

H ^ ■ A^ + 

^ 8192ai4 ^ 

/ m . 9 3a^m — 5/71^ . 4 5a^m — 28o?m^ + 27m^ . « 

+e ^A^ H A^ -7^ A^ 

' 2a2 16a6 64aio 

H ^ A^ + 

2048ai4 

2 /'a2 + m2 2 - 21m^ 4 , 50^ - 14a^m2 - 123a2m^ + 220m6 



-21a^ + 738a6m2 + 3080a^m^ - 19266a2m^ + 18445m^ .0 \ , 

+ A* + • • • (3.26) 

8192ai6 J 

In order to obtain the free energy, we consider the relation 

A J^F = 5.\/^A^ ^Aj + ^) = QsV^aCEXi) + C2 . (3.27) 

The vev A/) can be evaluated by looking at large z behavior of the planar resolvent 
Wofi{z) « gsVpN/z + gsVP{J2i ^i)/^'^ + ' ' ' in the loop equation (2.13). 

gs^/^A{Y,XI) = C2 + {^il- + (3.28) 
For /ii = /i3 = m, we have 

A— F = 2c2 = 2{A - m^) + A^ (3.29) 
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up to gl- By integrating in A, we obtain the free energy F = F^^ + eFo_i + e^Fo^2 + • • • for 
Np = 2 model. The planar part is given by 

A A 

Fo,o = 2{a' - m^) log A + A' + ^^^^ 

5a'^m? — 14a^m^ + 9m^ , r 

H A 

192aio 

5a^ - 252a^m2 + 1638a^m^ - 2860a2m6 + 1469^^ § 
+ 32768ai4 + ' ' ' ' ^^■^^> 

the half genus part by 

i^o 1 = ^A^ + s A^ A^ 

AS + ---, (3.31) 



8192ai4 

and the genus one part by 



0,2 



_ z + — , 

8a4 128a8 768ai2 



-2I0S + 738a^m'^ + 30800^^771^ - 19266a'^m'^ + 18445771^ 
32768ai 



+ ^ + • • • • (3-32) 



The Fo,o and Fq^i are the same as the known results [14] and [15], and the genus one 
correction Fq 2 exactly matches the corresponding part (B.16) of the genus one correction 
computed from the Nekrasov partition function. 

4 A^^ = 3 Model 

In this section, we consider the matrix model, which corresponds to SU{2) gauge theory 
with Np = 3. The action of the matrix model (2.4) is given by 

V{z) = Vo + eVi = 2^3log2; + 2777ilog(z - 1) - - + elogz. (4.1) 

z 



The function f{z) reads 
where ci , C2 and C3 are 



f{z) = — + + — , 4.2 

z z — 1 z-^ 



C3 = -5.V^E(a7)- ^^-^^ 



7=1 
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The condition V'{Xi)) = yields a constraint 

ci + C2 = . (4.4) 

From the large z limit of the planar loop equation in (2.13) with Wo,o{z) ~ gs^/JiN/ z and 
f{z) ~ (c2 + ca)/^;^, we have another constraint 

C2 + C3 = - (/is + mi)^ . (4.5) 

Since only one parameter among q's is independent, we choose C3 to parameterize the 
spectral curve. 

The planar spectral curve for Np = 3 model is of the form 

where the polynomial P4{z) is 

^4(2;) = 4m^/ - 4(5 - miA -m\ + ml^)z^ 

+{AB - 4miA + A2 - 4//3A)z2 + 2A(2/X3 - A)z + A^ , (4.7) 

with B = C3 — /isA + /z|. Henceforth, we take /is = m and ?ni = ?7t,oo = for simplicity. 
Then, P^iz) reduces to 

p^(z) = -ABz^ + (45 + A^ - AmK)z^ + 2A(2m - K)z + A^ . (4.8) 

The planar Coulomb branch parameter ag and the half genus part ao,i are 

2^/21^1^^"' ('-^^ 

with the derivatives of oq, 

9ao 1 / / A\ (iz daQ ^ f f dz 



dm 27ri J \ z J ^JPi ' dB l-ni J \ \fPi 
We will need the second derivatives of ao to compute the genus one part 00,2- 



(4.11) 
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Now, we will find the expression of ao,2 in terms of the derivatives of oq. The first term 
in (2.21) is given by 

i ^dz = i — \{z^ - z^)^j:r - 4z(3z - 2)^^ + = dz . 

2-niJ Syl 2-Ki J IQ V ^ p^l^ ^ ' p^l^ z-1 y/P^ j 

(4.15) 

The first integral on the right hand side of the above equation can be simplified using the 
fact 

S^ = -l ^ + f_L^ (4 16) 

p5/2 3dz2 /7^ + 3 3/2' ^^--^O^ 

and integrating by parts. 

The numerator of the last term satisfies the following relation. 
{z^ - z^)P'l = 2{z^ - z)Pi 

+2(45 + - 4mA)(z2 _ ^3) _^ 4A{2m - A){z - z^) (4.18) 

Then, we obtain the integral (4.17) after integration by parts 
1/1.3 Pa , 1 / 1 . . dz 

via (4.13) and (4.14). Plugging the above expression in (4.15), one has 



2^i 7 Sy^"^^ 2TTi J \12 y/Pi ^ A{z - 1) y/Pi) '^^ ^ 4 dB 

We give the computation of the integrals in the above equation in Appendix A. With (A. 7) 
and (A. 10), we obtain the first term of ao,2 as 

_}_ I yo_^ _ _}_ 1 gap m dap 
27Ti J 8y^ ^ ~ 12s "° 12 dB 12B dm 

- 4m-A — ^ + — - -125 + 2mA- A^' " 



8A' ' dm'^ 12A' ' dmdB 

-l(4i. + A2-4mA)^. 

(4.21) 
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The second term of ao,2 in (2.21) takes the form 



^ V u u 1^ ^ . A ^^ci _ ^z^ 2{z - z^) - (1 - z)) 



2m J 32y3 27ri J 4 ' ^ ' ^ 

A^a^gp A O^ao 1 a^ap 

8 4 5m5B 8 dra^ ' ^ ' ' 

The third term of ap^2 is given by 

1 / + 2V{' 1 / 1 - 2 , 



27ri / 8yp 27ri / 4Vi^ 

1 1 5ao m 5ap 

= -45"°+4 55+4B5^' (^-'^^ 

where we have used the result of (A. 7). Collecting all three terms (4.21), (4.22), and (4.23) 
of Op 2 , we find that 



ao,2 = -TTT^^o + o 



1 Oap m 9ap 



+ (_2m + A)— I + -(-65 + mA - 2A2) 



4A' 'am2 6A' ' dmdB 



+7(5-mA + A^)S. (4.24) 



The planar and the half genus contribution of the Coulomb branch parameter have 
been computed, which we reproduce here. [14], [15] 

A— m . 5 + 3m2 „ mihrr? + 5) . o 

^° = -^(^ + 45^-^45^^+ 25653 ^ 

A^ + -(^^^^ — " ' —A- 
6553655 

A6 + 0(A^)^ . (4.25) 



352 + 30m25 + 175m^ ^4 m(952 + 70m25 + 441m^; ^ 5 
163845^ ^ 6553655 

553 + 105^252 + 735m'^B + 4851m' 



104857656 

1 . m .9 5 + 3m2 . o 95m + 25m3 , . 

o-o 1 = ^A H 577rA ftttA H 7777^ — A* 

sVB 6453/2 51255/2 819257/2 

952 + 905m2 + 245m^ , , m (7552 ^ 4905^2 ^ i323m4\ 

A H ^A + • • • 

6553659/2 524288511/2 

(4.26) 
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The genus one part (4.24) is given by 



1 A — 51m^ 2 95m^ — 12 Bm ^ 

" ~ 1271 ^ 12B^ ^ 76855/2 + 1536^7/2 

1115^ + 1350W - 11725m^ ^^4 _ 7 (l55^m + 185Bm^ - WSSin^) 
196608^9/2 196608^11/2 

1315^3 + 32235^2^2 + 3814655771"^ + 635481r?i6 . « , , 

+ 25165824513/. + " " " " (^-^T) 



Inverting the vev a to get 5, we obtain 

„ , 9 mA + .9 — 677720^ + 570,^ . a 5a^m? — l^a?m^ + 9777^ . « 

B = [a^ \ ^A"' H ^ A* H -7^ A*" + 

2 32a2 8192a6 262144aio 

A 717 . 9 3a^m — . 4 ba^m — 2Sa?m? + 27777^ , 

+e A H A A + • ■ 

^ ' 4 32a2 ^ 4096a6 262144aio ^ 

9 / 1 771A — 7772 „ 5777^ — 3a 2777 . o 2a^ — 2ha?Tm? + 21777^ . . 

+e 1 1 A H A'^ A 

6 8a2 64a4 256a6 8192a8 

35a'^77i - 21002777^ + 207^7^ . . 

H T7^ A 



215a6 + 57190^^7772 + 7547702777^ + 6919337776 g , 



65536aio 
172 

4194304ai2 

The free energy can be computed from the derivative of F with respect to A. 



The planar and the half genus part of the free energy are integrated to 



Fo,o = {a" - m") log A + — + -^^A' 



H s ^ H T7^ ^ H ■ (4.30) 

32768 1572864 aio ^ ^ 



A ^ ,2 3a2777 — 5771^ . 5a^777 — 28a27n3 + 27771^ g 
= " 4 " 64^"^ + 16384 a6 ^ 1572864 aio ^ ^^'^^^ 
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The genus one part is given by 



i^0,2 = - - log A + — A + . A^ + A^ ^7,^:::^^ A 



,^4 

6'"*"'" ' 8a2'' ' I28a4 ' 768a6 32768a8 



H A 

327680aio 

215a6 + 5719a^m2 + 754770^^^ + 691933m6 
+ 25165824^12 ^ + ' ' ' " ^^■^^> 

Again, the free energy of the matrix model agrees with that of the gauge theory computed 
from the Nekrasov partition function in Appendix B.2. 

5 = 4 Model 

In this section, we consider the matrix model with action (2.5), which is related to the 
superconformal case of SU{2) gauge theory with four flavors. The Penner type potential 
reads 

Vo{z) = 2mo log z + 2mi log{z — 1) + 2m2 log(z — q) , 

yi(z) = iogz. (5.1) 

We evaluate the function f{z) in (2.12) as 

f{z) = — + - + , (5.2) 

z z — 1 z — q 

where the Cj's are given by 

/2mo + e\ /TTV^ / 2m/ 



Co = -gs 



Xi-l 



C2 = -..x/^E(^>- M 

From the equation of motion {J2j ^'(-^z)) = 0, it follows that 

2 

^Q = 0. (5.4) 

i=0 

And from the asymtotic behavior for 2: — )• 00 of the leading order loop equation in (2.13), 
one finds that the parameters satisfy 

ci+qc2=ml^- l^m-^ . (5.5) 

Therefore, we have a single parameter left to describe the spectral curve, which we take to 
be cq. 
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The planar spectral curve is 



(5.6) 



z'^{z — 1)^(2; — q)^ ' 

with Piiz) a polynomial of degree four. We will consider the case where the mass fii 
of all four hypermutiplets is equal to m such that the mass parameters rrii are set to 
mo = moo = and mi = m2 = m. Then, the polynomial P4:{z) becomes 

P^iz) = Cz^ + ((1 - qfm^ - C{1 + q))z^ + Cqz , 



(5.7) 



with C = qcQ. 

In this case, the Coulomb branch parameters up to half genus are 



ao 



ao,i 



1 



1 



liri J z{z — l){z 



Pa dz . 



(5.8) 



27ri 



m(l + q) 1 m(l — q)' 



dz 



{z-l){z-q)^^ 



dao Idao 
= — (1 + c?)^--^, 

where we have utilized in the last equality 

daQ 1 f 1 dz 



dC 2TTi 7 2 ^/R 
We compute the second derivatives of oq as 



dz 



{z-l){z 



'Pa 



dCdm 2-Ki 

d^ap _ J_ 
dm"^ 2TTi 



-m(l — q) 



r, dz 
z^ 



P 



3/2 



+ 



{l-qfz 



—m^{l — qpz'^ 1 
{z-l){z-q) '^/2 ' 



dz , 



(5.9) 

(5.10) 
(5.11) 

(5.12) 
(5.13) 



Now, we calculate the genus one part ao,2 in terms of the derivatives of ap. The first 
term of 09,2 in (2.21) becomes 



Im J 



dz 



1 



z{z-l){z-q) P; 



3z2 



2{l +q)z + q P^ 



32 



P 



5/2 



p; 



,3/2 



(3^2 -2(l + g)z + g)2 1 



dz . 



Sz{z-l){z-q) /f^y"'~' ^^'^^^ 

Using (4.16) and integrating by parts, the first integral on the right hand side can be 
written as 



1 



^z-{l + q) 1 ^z{z-l){z-q) P^ 



12 



'Pa 



48 



P 



3/2 



dz . 



(5.15) 



- 15 - 



This can be written using (A. 14) as 



I r(3 z m^{l-qfq z \ 1 dao m'^{l-qfdao 



2vri J yS^/P^ 8 p3/2 J 4^ Qc 6C dC 

m{m^{l-qf -C{l+q)) d'^ap 

6C dCdm ' 

The second term in (5.14) is integrated by parts to 



(5.16) 



1 // 3 1,, ^\ dz 1 / — 3 2; , ,9ao 

^' z + -{l + q)] ^= = — d) —dz + {l+q)-^. (5.17) 



2-Ki J \ 2 2' 7 2m J 2 y/P^ ' ' dC 

The third term in (5.14) can be written as 



27riJ8\ ' z ^ ^' {z-l){z-q)J ^ 

1 /9 z , 1 f m'^n- qf -C(l + q) qz , 

zdz - (B TTTTCtZ 



27ri / 8^/ft 27ri / 8 p3/2 



4 ^TTi J 6 

_3, s^^Qq gg 19^00 1 dap 

V '^^dC 2m{l-q)^dCdm 8mdm' ^' ' 

where (A. 21) has been used. Adding (5.16), (5.17), and (5.18), the first term of 09,2 in 
(5.14) becomes 

,,'2 



2vri J 8yl 27ri J 8 p3/2 



4 

?n2(l-g)2 5ao m(m2(l -g)^ -C(l + g)) d'^ap 



6C dc 6C acam 

Cq d'^ap 1 dap 

+ ^izr- (5-19) 



2m(l — q)"^ dCdm 8m dm 
For the integration of z/P^ in the above equation, we use (A. 17), which results in 

1 I y'p r 1/-, , ^^ap 1+q , 2/-, ^2 , "9^00 

m"^ {1 — q)'^ dap m(m^(l — g)^ — C(l + g)) d'^ap 



6C dC 6C dCdm 

I Cl d'^ap I 1 dap 
2m{l — qY dCdm 8m dm 
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Next, we compute the second term of ao,2 in (2.21). 



27rj / 32y3 

,2 / (^ _ „^4,3 



/ 



1 (9^ao 1 (?ao 



8 9m-^ 8m dm 



1 9^ao 1 dao m^(l + g)^9ao 
~8 9m2 8^ 5m 4C 9C 



The third term of the genus one part of the Coulomb branch parameter (2.21) is 
calculated to 

1 / Vf + 2V{' , 
'* -dz 



27ri / 8yo 

1 f —1 / ,^ , Q\ dz 



27ri / 8 V ' ' zJ 

?Ti2(l - qf - C7(l + q) dao ^ {^^{1 - q? - C{1 + q)f d^ao 



8C dC 8Cm{l-qY dCdm 

Cq d'^ao _^ l^--^ _^ Cq d'^ao 



2m(l - g)2 dCdm 4 ' dC 2m(l - q)'^ dCdm 

+ 2^/l("^'(^-^)'-^^^ + ^))^^^ 
m^l - qf dao ^ (^^(1 - q? - C{1 + q)f d'^ao 



4C dC 4Cm(l-g)2 SC^m 



m(l — 5)2 dCdm 



(5.22) 



where (A. 20), (5.10), and (A. 21) have been used in the second equality and (A. 17) was 
used in the last equality. From the results of the three terms of 00,2 in (5.20), (5.21), and 
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(5.22), we find the genus one part of tlie Coulomb branch parameter (2.21) 
1 .5ao , m'^{l-qfdao Cq d'^ao 



dC VIC dC 2m{l-qYdCdm 

, ^ l + g I 2/1 n2 n^^ 1 ^^ ^^"0 ^ ( 2/. n2 n^^ , \\ ^^"0 
+ ^7^ IT? ("^ (1 - ^) - C^il + g)) a^n^ - 7^("^ (1 - - C"!! + 9))- 



8?Ti (1 — (/)^ dCdm 6C dCdm 



8 5m2 4C 9C AC{l-q)^^ ' ' ^" dCdm 

_^m(l+_g)^_a^ _2^g(l + g) S^oq 



4 (1 - g)2 5Cam (1 - qY dCdm 

I (m2(l-g)2-C(l + g))2 d'a, 

4Cm(l - g)2 5C5m ' ^ ^ 



We give a few leading order terms in m of the planar and the half genus contribution 
of the Coulomb branch parameter. [14], [15] 



ao = ivC /io(g) - hi{q)— 5—7^ ^777 ^771 + C 



C 3 5 C3 7 C4 V C5 

771 777/^ 777^ 777*^ \ 

ao,i = i l^giiq)-^ + 92{q)-^ + 53(g)^j;^ + + ■■■ j , (5.24) 

where the functions hj(g)'s are obtained from the expansion of the hypergeometric function, 

= i + i" + ik"' + 512"' + siss"* + "^""'^ '^'^''> 
3 27 27 3 27 

'•^<") = 8 + 32" + 512"' + 2048"' + I3IO72"' + °<"''' <=^^^^' 
, , , 5 125 1125 2 125 , 125 . 

"'<"* = 16 + liT" + I024" + 4096" + 262144" + °(" >• '^'^''^ 
, , , 35 1715 42875 n 42875 o 42875 4 

= 128 + + + 3^8^ + 5097l5^'^ + ^ ^'-''^ 



and the function ^'^(g) is given in terms of /ij(g)'s by 



2/ii(g)-(l + g)/io(g) 
91(g) = 7; > (5-30) 



92ig) = (5.31) 

/ N 6/t3(g) - 5(1 + g)/t2(g) , . 

ff3(g) = 77; , (5.32) 



9A{q) = ^' (5.33) 





2 


4/12(9) 


-3(l + g)/ii(g) 




6 


6/i3(g) 


-5(l + g)/i2(g) 




10 


8/14(9) 


-7{1 +q)hs{q) 
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The genus one part of the vev (5.23) is computed to 

(1 TTi^ Tfl^ TUp \ 

^1(9)-^ + f^2iq)-^ + hiq)-^ + ki{q)-^ + ■■■ j , (5.34) 

where ki{q) is related to /ii((?)'s as 

ki{q) = -^{-{q + ^)ho{q) + 2hi{q)) , (5.35) 
lb 

k2{q) = W + 4? + ^)ho{q) - I7{q + l)/ii(g) + 12/i2(<?)) , (5.36) 

ksiq) = (12((?2 + 4g + l)hi{q) - 3l{q + l)/i2((?) + I8hs{q)) , (5.37) 

A:4(g) = (20((?2 + 4g + l)/i2((?) - 45((? + l)/i3((?) + 24/i4(g)) . (5.38) 

With the exphcit resuhs of the Coulomb branch parameter, we solve the equation a = 
ao + eao,i + e^ao,2 for C. 

^ _ z, _J: 2/ii(g) 2ho{q)h2{q) - Shijqf 

" h^{q? ho{q) + 3 

Who{q)hi{qf - Who{qfhi{q)h2{q)+2ho{qy^h3iq) 



5 a6 

^^^^ ' 2(l + g)/io(g) -4/ii(g) ^ 8/io(g)/i2(g) - 12/ii(g)^ 
/lo(q') 3 

mo(g){5/ii(g)3 - 5/io(g)/ti(g)/i2(g) + /to(g)'/t3(g)} 

5 



+ 



+ ••• 



+ 6^ 



(l + g)/io(g)-2/ii(g) 
8/io(g) 

15/ti(g)' + /io(g)'(-l + qf - ho{q){6h2iq) + 5/ii(g)(l + g)) 

12 a2 

{81/11 (g)^ - 3/io(g)/ii((?) (23/12 ((?) + 5(1 + q)hiiq)) 



12 

+/io(g)2(9/i3(g) + 10(1 + q)h2{q))}'^ + ■■■ 



(5.39) 



Now, we can compute the free energy by 



|--».v/5(Ex^)-^- (-0) 



From (5.5), C2 is related to C = qco, thus we obtain 

±F = - 
dq 1 



l-F = ( 4m2 - . (5.41) 
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The Fo,o and Fo,i are given by 
-Fb,o = (a - m ) log g H 



+ 



2a2 

(I3a^ + lOOa^m^ + 22a'^m^ - I2a'^m^ + 5m^) 2 



192ai 



+ TT^^o ^ 



+—^—rj (2701a''' + 26440a"''m^ + 7164a"^m'* - 9000a'^m° 
327680^"* 

+12190a^m^ - 13384a^m^° + 10908a^m^2 - 5720a^m^^ + 1469m^^) + ■ ■ ■ , 

(5.42) 

^ 2m(a2 + m^) 9a^m + lla^m^ — 90^/71^ + 5m'' n 

-Pb 1 = ^ o -Q a Q 

2a? 16a6 ^ 



96aio 

1 



1257a^^m + 1791a^2^3 _ 3375^10^5 _^ 6095a^m^ - 8365a^m^ 



4096ai4 

+8181o''m^^ - bOOba'^m^^ + 1469m^^) + • • • . (5.43) 

Finally, we obtain the genus one free energy i<b,2- 

+ Ga^rr? + m"^ 9a^ + 640*^7712 - 70a^m^ + AOa^m^ + 21m^ 2 
- 8^^ ^ + 128^8 

190^2 _^ i47aio^2 _ 3ooa8m^ + 4700^777^ - 357a''7778 + 390^777^° + llOm^'^ 3 
^ 384^12 

H 3— (I257a^^ + 102760^^777^ - 287760^^777^ + 67460a^°77i^ - 105630a*^777^ 

32768ai6 ^ 

+103308a^777^° - 43120a^777^^ - 150280^777^'' + 18445777^'^) + • • • . 

(5.44) 

The genus one part is identical to (B.6), which is calculated in A/" = 2 SU{2) superconformal 
gauge theory with four flavors. 
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6 Summary and Discussion 

In this paper, we have investigated the /^-deformed matrix models with logarithmic po- 
tentials, which have been suggested to explain the AGT relation. Specifically, we have 
computed the genus one part of the free energy in the matrix models, which describe 
Af = 2 SU{2) gauge theories with Np = 2, 3, and 4 flavors. We have checked that the 
results obtained from the matrix model nicely match the free energies computed from the 
Nekrasov partition function. 

It would be interesting to generalize the duality to SU{Nc) gauge group, in which 
case a multi-matrix model should be considered. [16] [17] One can also consider extending 
the results to higher genus parts of the free energy. [33] , [34] , [35] The /3-deformed matrix 
model is also related to the refinement of the topological B-model. [36] The study of matrix 
model along this line would be highly interesting. 
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A Calculation of Integrals 

In this appendix, we show the details of calculations of the results used throughout the 
paper. 

A.l Integrals for Np = 3 Model 

Consider the oq in (4.9) in the form of 

2-Ki T 2\z-l z z- 



We will rewrite the three terms in the above equation in turn. Through the relation 

-l- = {-iBz'' + {-AmA + +A^)z-A^)^, (A.2) 
z — 1 

the first integral is written as 

2 dz 



— i -^Jpldz = — li-ABz'^ + {-AmA + h?)z- k , ^ 
27ri J z - 1^ 2-Ki J ^ ^ ' ' 
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where the following condition has been used to eliminate term. 



1 

2^ 



dz 



P' 

dz 



2vri J 2/^ 



— (fi-QBz'^ + UB - 4mA + A^)^ + A(2m - A)) 
2m J 



dz 

7k 



(A.4) 



Similar manipulations give the second term in (A.l). 



im J z 



1 

27ri 

1 

2^ 



-4Sz^ + (45 - 4mA + k^)z + 2A(2m - A) + — 

z 

-(45 - 4mA + A^) + -A(2m - A) + — ^= 

3 3 Z J y i-*4 



A2\ dz 



(A.5) 



The third term of (A.l) can be written as 



1 

2TTi 

1 

27ri 



P' 



z \ 2JPa 



-QBz + (45 - 4mA + A^) + ^t^^!IlJzR\ ^ . (a.6) 



'Pa 



Adding (A. 3), (A.5), and (A.6), one finds (A.l) to be 



(A.7) 



where (4.11) has been used. 



Next, we rewrite the following integral 



1 dz 
2m J z-1 JVa 



— ([{-ABz'^ + (-4mA + A'^)z - A^ 
m J 



+ — ([{-{4:B + 4?nA - A'^)z - A^) 
2m J 



2m 

2B 52 ao 
A dmdB 

^^ + ^(45 + 4mA-A^)^ 
A dmdB + 2A2 + ^""^^ ^ 5m2 

H : * (-45 - 4mA) ^ 



dz 



dz 



P 



,3/2 



2m 



p: 



,3/2 ' 



(A. 



where (4.13) and (4.12) have been used in the second and third equality. The integral in 
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the last line can be obtained from the relation 



1 f P' 

= — (t -^dz 



2TTi T p3/2 

dz 



— (p (-UBz^ + 2UB - 4mA + A^)^ + 4mA - 2A' 
2vr^ / 

6B d'^an 1 , „ . . .9s5^ao 



p3/2 



Thus, we find the integral (A. 8) to be 

^ +ttt^(-4"^A + A2)— ^. (A.IO) 



27ri / z - 1 ^/A A dmdB 2h? ^ ' dm? ' 



A. 2 Integrals for Np = 4: Model 

First, we consider the integral in (5.15) 



1 f 1 



27ri 



j^{P,-m\l-qfz')^^dz 



1 / /6Cz-2C(l + g) + 2m2(l-g)2 m2(l - g)2 ^2p^/ 



dz . 

(A.ll) 



We rewrite the second term in the last line via 

z^P^ = 6P4 - 4(m2(l - qf - C{1 + q))z^ - 6Cqz , (A.12) 
and (5.12) to obtain 

1 f -m^(l-qf z'^Pi' 1 f f -6m^(l - q)^ 1 .2 ^ \ , 

^ = TT- <^ 7^ 7^+6m2(l-g)2g— - 



27ri/ C p3/2 27ri J I C '^^'^p3/2 
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Substituting the above equation into (A. 11), we get 



If/ , , 4m2(l-fl)2\ dz 1 / x9 z , 



I f [ 6z ^ n9 \ , Z^./-, N 8m^(l — dao 



27ri/ ' ^' p3/2 y ' ' (J J Qfj 

where (5.10) has been used in the last equahty. In the above expression, we need to compute 
the integrals of z/ and z jP^ , which we do in the following. 

So, we rewrite da^/dC of (5.10) to obtain the integral of z/Pj through 

zP'^ = 3P4 - ("1^(1 - q? - C(l + g))z2 - 2Cqz , (A.15) 

such that 

dan 1 f dz 1 f z P'. ^ 

-dz 



dC 27ri J 27ri J 4 p3/2 



_L ^ ^ (3P4 - - qf - C{1 + q))z^ - 2Cqz) 



p3/2 



^dao m\l-qf -C{l + q) d^ap 1 f -Cqz dz 

2dC 2m(l - g)2 dCdm 2m J 2 pV^ ' ^ ' ' 



Therefore, we obtain 



27ri / p3/2 Qc m{l - qf dCdm 

Next, we will compute the integral of zj^fP^ using the relation 

+ ^"^'^^ - '^^^ - "^^^ + + (n^^d - qf - en + ,)),. . (A.18) 
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It follows that 

1 f z , I r z'^ Pi 

zdz — '* 



2-Ki J ^fPi 2m J 4 p3/2 



1 r ( 3z m2(l-g)2-C(l + g) 1 Cqz^ 



2TTi J \A^/P^ 4:C ^ 2Pf/^ 

(m^(l-g)^-C(l + g))^ z^ m\l - - C{1 + q) qz \ 

1 rS z m^{l-q)^ -C{l + q)dao ^ Cq d'^ao 



2TTi J 4 /f^ 2C dC m{l - qf dCdm 

im\l-qf-Cil+q)f d^ap 
2Cm(l - g)2 dCdm 

m^{l-qf -C{l + q) / dap m^{l - qf - C{1 + q) d'^ap \ 
^ AC \dC^ m(l -g)2 dCdm) ' 

(A.19) 

where (A. 17) has been used in the last line. Thus, we obtain 
1 / ^ dz- "^^C^-l)^ -C{'^ + q)dap 



2m J y/Pi C dC 

(^(l_g)2_c(i + g))2 Q2^^ ^ ^cq d^ap 
Cm{l — q)^ dCdm m(l — g)^ dCdm ' 
The integral used in (5.18) is computed as 



— f{-2Cz^-im\l-qr-C{l + q))z) 
■J P4 

AC d'^ap 1 f , 2/1 \2 \N dz 



^ f{m'{l - qf - C{1 + q))z^^ , (A.21) 

^ Pa 



m(l — g)2 dCdm 2-Ki J 



where in the second equality the relation 

1 f P' 

^=2mfjt2^' 



was used to replace the constant term. 
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B Nekrasov Partition Function 



We give a brief summary of the Nekrasov instanton partition function and present relevant 
results obtained from them. [37], [38] 



B.l Nf = 4: theory 

The instanton part of the Nekrasov partition function for M = 2 U{2) gauge theory with 
four anti-fundamental hypermultiplets is of the form 

4 

^inst = ^ q^^^Z^cc{a, Y) JJ ^afund(a, Y, ^J-i) , (B.l) 

where Y = (11,12) is a pair of Young diagrams, a = (01,02) a pair of Coulomb branch 
parameters and /Xj denotes the mass of the hypermutiplet. The vector mutiplet and anti- 
fundamental hypermultiplet contributions are given by 

2 

Zvec(a,r)= n lliav-^iLY,is)+e2iAYM + ^)y' lliaji+ei{LY^it) + l)-e2AYMy\ 

i,j=l s£Yi t£Yj 



2 

^afund(a, ^, /^) = n n + ei("i - 1) + ^2in - 1) + /i) , (B.2) 

i=l seYi 

where Uij = Ui—aj, and the leg- length -Lyj(s) = A^ — m and the arm-length ^y-(s) = Xm — n 
are defined for a box s at position (m,n) in a Young diagram Yi = (Ai > A2 > ■ ■ ■), with 
its transpose Y^ = (A'^ > A2 > • • • ). 

To compare with the results from matrix models, we consider SU{2) gauge group with 
the Coulomb branch parameter a = (a, —a), noting that the contributions coming from the 
C/(l) factors are irrelevant for genus one correction. We define the free energy as 

J" = -6162 log Zinst 

= Y.gfe'Tk,i.. (B.3) 

k,l=0 

For the case of fii = m, which we considered in the matrix model computation, we calculate 
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the free energy up to gj order. 

+ 6a^m^ + I3a^ + lOOa^m^ + 22a^m^ - 12a^m^ + 5m^ 2 
- 2^^ + 64^6 9 

23a^2 _^ 204ai°m2 + Sla^m'' - 48a''m'' + Aba^m^ - 28a'^m^^ + 9m^^ 3 

192^10 ^ 

32768ai4 



+ 



+ ^4 (2701a^'' + 26440o^^m2 + 7164a^^m'* - 9000a 



+ 12190a^m'^ - 13384a^m^° + lOOOSa^m^^ - 5720a'^m^'^ + 1469771^*^) 

(B.4) 



= 16^^ 

38ai°m + Sla^m^ - 72a^7n^ + 90a'^m'^ - TOa^m^ + 27m^'^ 

1 



-q 



3 



{l257a^Sn + ITgia^^n^ - 3375a^°m^ + 6095a^m^ - 8365a^m^ 



40960^4 

+8181a^m^^ - SOOSa^jn^^ + 1469m^^) q'^ + O(g^) , (B.5) 



+ ea^m^ + 9a^ + 640*^771^ - 70a^m^ + AOa^m'^ + 21m^ 2 
- 8^^ + 128^^ ^ 

19ai2 ^ I47ai0m2 - 300a8m'' + A70a^m^ - 357a^m^ + 39a^m^° + llOm^^ 



^ 384ai2 



+ ^--r(1257ai6 + 10276a^^m2 - 28776a^2m^ + 67460a^^m^ - 105630a^m^ 

32768ai6 ^ 

+103308a^m^° - 43120a^m^2 - 15028a'^m^'^ + 18445?n^^)/ + 0{q^) , (B.6) 

m2(— a^ + m^)^ 2 m2(— + m2)^(3a^ — 3a2m2 + 8m^) 3 
= 32^^ + 96^^^ ^ 

m2(-a2 + m2)3(235aS - 432a^m^ + 1486a'^m^ - lesea^m^ + \U7m^) 4 5 



-5^ + 0((f 



8192ai6 

(B.7) 
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B.2 Nf = 3 theory 



The instanton partition function for this case is 

3 

Z-mst = ^ AIj^' Zvcc(a, Y) Y[ ^afund(a, Y, fii) , (B.8) 



Y 



1=1 



where A3 is the dynamical scale of the gauge theory. As in section 4, we set /xi = /U2 = 
and /U3 = m. The free energy is computed as 

= y ^ + ^4^^ + 32768ae ^ + " " " ' (^"^^ 

^o, = -A-^A^ + "^('"'"'f^ A^ + ---, (B.IO) 
4 64a2 16384 a6 ' ^ ^ 

m '^^~'^^a2 m(— Sa^ + Sm^) 3 
-^o,2 = ^A + ^^A + ^^^-^ A 

= -^28^^ + 32768^^ ^ + ' ' ' ' ^^'^^^ 



B.3 Nf = 2 theory 

In this case, the instanton partition function is of the form 

Z\nst = ^ Ag'^'Zvccla, Y) ^afund(a, /^i) • (B.13) 

Y 1=1,3 

By setting /^i = ^3 = m as in section 3, we find the free energy as 

^ a? + m? .9 — Qa?m? + 5m'* . a ^a^m? — lAa^m^ + 9m^ . e 

Jj) = n — H ^ A* H -7^ A'' 

2a2 64a6 192 a^o 

5a^ - 2b2a^m? + IGSSa^m^ - 28600^7716 + 1469771^ . § ^. 
+ 32768^^^ ^ + • • • ' ^^-^^^ 
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Tq 1 = ^A^ + — 7 -h^ T?^ ^A*^ 



^TTTTT^i ^A8 + ..., (B.15) 



m 

^~ 8192 



+ „ - 21m^ . . 5a^ - 14a^m2 - I22,a^m'^ + 220771^ , « 

JJ) 2 = 3 — A"' 3 — A* H A'' 

8a4 128a8 768oi2 

H A^ H , (B.16) 

32768ai6 ^ > V ^ 



- 3a^m^ + 2m^ 4 lla^ - 27 a^m'^ + 16m^ q 
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